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The  Dopant  Density  and  Temperature  Dependence 

of 

Electron  Mobility  and  Resistivity 
in  N-Ty-pe  Silicon 

by 

Sheng  S.  Li 


Abstract:     Traditional  analysis  of  electron  mobility  in  n-type  silicon  ne- 
glects the  effects  of  electron-electron  scattering  and  scattering  anisotropy  in 
the  mobility  calculations.     As  a  result,  theory  fails  to  conform  with  experiment 
when  dopant  density  exceeds  2  x  10^^  cm~^. 

In  this  work,  an  improved  theoretical  model  for  computing  mobility  and  re- 
sistivity as  functions  of  dopant  density  and  temperature  has  been  formulated  for 
n-type  silicon.     The  model  has  been  applied  to  phosphorus-doped  silicon  for  dopant 
densities  from  10^^  to  10^^  cm~^,  and  temperatures  between  100  and  500  K. 

The  electron  mobility  was  calculated  analytically  by  appropriately  combining 
lattice,  ionized  impurity,  and  neutral  impurity  scattering  contributions.  The 
effect  of  electron-electron  scattering  on  both  lattice  and  ionized  impurity  scat- 
tering mobilities  was  incorporated  empirically  for  dopant  densities  greater  than 
2  X  10^^  cm~^.     Additionally,  the  anisotropic  scattering  effect  was  also  considered 
in  the  mobility  calculations. 

Theoretical  calculations  of  resistivity  and  ionized  donor  density  were  made 
for  phosphorus-doped  silicon  for  temperatures  between  100  and  500  K. 

Resistivity  measurements  on  seven  phosphorus-doped  silicon  slices  with  do- 
pant densities  from  1.2  x  lO-^'*  to  2.5  x  10^^  cm~^  were  carried  out  for  tempera- 
tures between  100  and  500  K.     Electron  mobility  at  300  K  was  deduced  from  resis- 
tivity and  junction  C-V  measurements  for  dopant  densities  from  10^'*  to  10^^  cm~^. 

Agreement  between  theory  and  experiment  for  both  electron  mobility  and  re- 
sistivity of  phosphorus-doped  silicon  was  within  ±7  percent  in  the  range  of  do- 
pant densities  and  temperatures  studied. 

Key  Words:     Dopant  density;  electron  mobility;  ionized  impurity  scattering 
mobility;  lattice  mobility;  neutral  impurity  scattering  mobility;  n-type  sili- 
con; resistivity;  scattering  mechanisms;  temperature. 


1 .  INTRODUCTION 

Although  a  considerable  amount  of  work  has  been  published  which  deals  with  electrical  prop- 
erties and  mobilities  for  silicon,  it  is  evident  that  traditional  analysis  of  electron  mobil- 
ity in  n-type  silicon  fails  to  predict  correctly  the  mobility  values  when  dopant  densities 
exceed  2  x  10-^^  cm~^.     The  discrepancy  may  be  attributed  to  the  inadequacy  of  the  existing 
theoretical  models  in  the  high  dopant  density  range.     For  example,  the  effects  of  electron- 
electron  collisions  and  scattering  anisotropy  on  both  lattice  and  ionized  impurity  scatter- 
ing mobilities  have  been  neglected  in  the  traditional  analysis. 

In  this  work,  an  improved  theoretical  model  for  computing  electron  mobility  and  resistivity 
as  functions  of  dopant  density  and  temperature  for  phosphorus -doped  silicon  is  described. 
The  calculations  cover  the  range  of  dopant  densities  from  10^^  to  10^^  cm"^,  and  tempera- 
tures between  100  and  500  K.     Detailed  comparison  of  the  existing  scattering  theories  and 
mobility  formulations  for  n-type  silicon  is  made  in  section  2  and  Appendix  A.  Calculations 
of  electron  mobility  and  resistivity  for  phosphorus-doped  silicon  are  described  in  section  3. 
Results  and  discussions  are  given  in  section  4.     Section  5  presents  the  conclusions  from  this 
work. 
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2.     SCATTERING  MECHANISMS  AND  MOBILITY  FORMULATIONS 


In  this  section,  a  general  review  of  the  existing  scattering  theories  and  mobility  formula- 
tions for  n-type  silicon  is  given.     Scattering  processes  due  to  lattice,  ionized  impurity, 
and  neutral  impurity  that  contribute  to  the  total  electron  mobility  in  n-type  silicon  are 
described.     The  mobility  formulation  for  each  of  these  scattering  processes  is  described. 
Quantitative  comparison  among  the  mobility  values  calculated  from  the  different  theoretical 
models  for  lattice,  ionized  impurity,  and  neutral  impurity  scattering  is  made.     The  effects 
of  electron-electron  collisions  and  scattering  anisotropy  on  both  lattice  and  ionized  impu- 
rity scattering  mobilities  are  discussed. 

2.1.     Lattice  Scattering 

Calculations  of  lattice  scattering  mobility  for  n-type  silicon  have  been  reported  by  Long 
[1],  Norton  et  dl.    [2],  and  others  [3,4].     All  these  calculations  are  based  on  the  general 
treatment  developed  earlier  by  Herring  and  Vogt  [5]  for  multivalley  semiconductors  such  as 
silicon  and  germanium.     The  lattice  scattering  mobility  for  n-type  silicon  is  mainly  due  to 
the  intravalley  acoustical  phonon  and  the  intervalley  phonon  scattering  (i.e.,  the  so-called 
f-type  and  g-type  phonon  scattering)   [2].     The  lattice  mobility  due  to  acoustical  phonon 
scattering  has  been  discussed  by  Brooks   [4],  using  the  concept  of  deformation  potential  in- 
troduced first  by  Shockley  and  Bardeen  [6],  for  semiconductors  with  spherical  energy  sur- 
faces.    The  reciprocal  relaxation  time  for  this  case  is  [4]. 

T^l  =  (8^3/h^)    (2m*) (kT/pu^2)   ^^^2  £1/2  (2.1) 

which  predicts  that  t  is  proportional  to  E~^/2_  xhe  lattice  mobility  derived  from  eq  (2.1) 
is*  [4] 

2(2^)1/2  qh"*  pu  2 

M      =  ^  <T     >  =  —  (2.2) 

La      m*      La        3^^5/2  ^^^^3/2  g_^^2 

Equation  (2.2)  predicts  that  p  varies  as  T~^*^  for  acoustical  phonon  scattering.     For  sili- 
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con,   (m*)~^'2  is  replaced  by  -r-  (m  m  2)-l/2  (m        +  2m  ~^)  and  (m  /m*)^/2  =  20.4;  pu  2  = 

.3  ^  t                   X/               t                    o  /v 

1.97  X  10-^2  dynes/cm2;  E^    =  12.80  eV  is  the  deformation  potential  constant  [4].  Lattice 
mobility  described  by  eq  T2.2)  assumes  that  all  the  scattering  events  are  due  to  dilation, 
and  that  scattering  between  different  ellipsoidal  energy  surfaces  (i.e.,  the  intervalley 
phonon  scattering)  has  been  neglected.     The  experimental  observations  of  electron  mobility 
for  intrinsic  silicon^  however,  shows  a  stronger  temperature  dependence  of  the  lattice  mobil- 
ity (i.e.,  Uj   ~  x~2.3-0.2)  than  that  predicted  by  eq  (2.2).     To  correct  this  discrepancy. 
Herring  [7]  has  also  considered  the  optical  mode  and  the  intervalley  phonon  scattering  for 
n-type  silicon.     Long  [1]  has  extended  the  expression  to  a  general  case  in  which  there  are 
several  intervalley  phonons  of  different  energies  and  scattering  strengths  and  in  which  the 
acoustic  phonon  scattering  may  be  anisotropic.     Herring  and  Vogt  [5]  have  shown  that  it  is 
legitimate  to  describe  the  scattering  by  a  relaxation  time  tensor  diagonal  in  the  principal 
axes  of  an  energy  spheroid.     There  are  two  relaxation  times  for  n-type  silicon,        and  t^. 
The  reciprocal  relaxation  time  component  for  this  case  is  thus  given  by  [1] : 


E     1/2    T  1  i 

o  010 


3/2 


-1 


(E/kT   .  +  1)1/2 
ci 


exp  (T  ./T)  -  1 
ci 


(2.3) 


(E/kT   .  -  1)1/2 
ci 


1  -  exp(-T^./T) 


*The  average  relaxation  time  <t>  for  the  nondegenerate  case  is  <t>  =  /°°  t  e3/2  exp (-E/kT)dE/ 
r  e3/2  exp(-E/kT)dE  where  x  is  given  by  (2.1).  ° 


The  W.     and  W.  measure  the  strength  of  coupling  of  the  electrons  to  acoustic  and  to  inter- 
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valley  phonons,  respectively,  and  the  subscript  i  ranges  over  all  the  intervalley  phonons  of 
different  energies.     The  subscript  a  is  to  be  replaced  by  "t"  or  "I"  depending  on  the  direc- 
tion relative  to  the  spheroid  axis.     The  temperatures  T,  T  . ,  and  T    are,  respectively,  the 
temperature  at  which  the  relaxation  time  is  to  be  determined,  the  cBaracteristic  temperature 
of  the  ith  phonon,  and  a  reference  temperature  chosen  to  fix  the  magnitude  of  x.     The  lat- 
tice scattering  mobility  for  n-type  silicon  can  be  calculated  from  eq  (2.3),  by  use  of  the 
following  expression: 

The  procedure  for  computing        from  eq  (2.4)  is  to  substitute  eq  (2.3)  into  eq  (2.4)  and  to 
adjust  the  intervalley  coupling  constants  W.   to  give  the  best  fit  of  the  resulting  expres- 
sion to  the  observed        versus  T  curve.     To  simplify  the  theoretical  calculations  of  U-^  by 
eqs  (2.3)  and  (2.4),  Long  [1]  has  assumed  that  the  possible  contributions  by  all  of  the  vari- 
ous intervalley  phonons  can  be  approximated  by  just  two  phonons  of  energies  equivalent  to 
temperatures  of  630  to  190  K.     This  approximation  is  reasonable,  since  all  of  the  intervalley 
phonons  lie  within  two  rather  narrow  energy  ranges,  and  these  phonons  approximately  represent 
average  energies  of  the  two  ranges.     The  scattering  anisotropy  due  to  the  anistropic  effec- 
tive masses  m^  and  m^^  is  handled  by  properly  choosing  the  ratio  of  coupling  constants  W.  and 

in  the  mobility  calculations.     Long  [1]  has  found  that  choosing  the  ratio  of  W^^/Wy^^^  =  2.0 
anS  ^2^^^t  ~  ^'^^  yields  the  best  overall  fit.     The  dashed  line  in  figure  1  shows  the  results 
of  his  calculations  for  temperatures  between  100  and  500  K.     Norton  et  at.    [2]  have  recently 
reported  a  similar  calculation  (but  assuming  four  intervalley  phonons  with  phonon  tempera- 
tures       =  540  K,  Tj  =  690  K,  T^  =  190  K,  and  T^  =  307  K)  for  phosphorus-doped  silicon. 
Their  results  are  also  shown  in  figure  1  by  the  solid  line.     It  is  noted  that  both  curves 
are  in  good  agreement   (within  ±3  percent)  for  temperatures  between  100  and  500  K.     For  a  more 
detailed  discussion  of  their  individual  treatments,  the  readers  are  referred  to  the  original 
'    papers  by  Long  [1]  and  Norton  et  dl.    [2],  respectively.     We  have  used  the  lattice  mobility 
data  of  Norton  et  at.   [2]  for  the  present  mobility  calculations. 

2.2.     Ionized  Impurity  Scattering 

For  lightly-doped  n-type  silicon  the  main  contribution  to  the  electron  mobility  comes  from 
lattice  scattering.     As  the  dopant  density  increases   (or  temperature  decreases),  the  role  of 
impurity  scattering  becomes  more  important.     Theories  for  ionized  impurity  scattering  in 

I    semiconductors  have  been  developed  by  Brooks  [4]  and  Herring  [7],  Dingle  [8],  Samoilovich  et 

\   at.    [9],  and  Luong  and  Shaw  [10]. 

The  Brooks-Herring  formula  is  based  on  the  important  assumptions  that  the  Born  approximation 
applies,  the  relaxation  time  is  scalar,  the  energy  surfaces  are  spherical,  electron-electron 
interactions  are  negligible,  and  impurity  cell  effects  can  be  ignored.     The  expression  for 
the  relaxation  time  derived  from  these  assumptions,  is  [7]: 

(2m*)l/2  e  2  e3/2 

T    =   X  10-6  (2.5) 

irq"*  N^  G(b) 

where  G(b)  =  In  (b  +  1)  -   [b/(b  +  1)],  (2.6) 

24™''  e  (kT)2 

b  =  — ^   X  10-6,  and  (2.7) 

q2h2n' 


n'  =  n  +  n  (1  -  n/N^^) ,  for        =  0  (2.8) 

Equation  (2.5)  predicts  that  t  is  proportional  to  E^^2.  xhe  ionized  impurity  scattering  mo- 
bility may  be  derived  from  eq   (2.5)  with  the  result  that  [7] 


(2.4) 


3 


Figure  1.     Lattice  scattering  mobility  as  a  function  of  temperature  for  n-type 
silicon  calculated  by  Long  [1]   (dashed  line)  and  Norton  et  at.    [2]   (solid  line). 
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^TCRH^  =  — ;  ir~>  10-2  (2.9) 

^^^"-^        ^3/2  q3m  1/2  N  G(b) 

Equation  (2.9)  takes  no  explicit  account  of  any  anisotropy  in  the  ion  scattering.  As  a  re- 
sult, it  usually  overestimates  the  mobility  values,  particularly  in  the  high  dopant  density 
range  where  impurity  scattering  is  dominant. 


The  difficulty  in  using  eq  (2.9)  for  computing  the  ionized  impurity  scattering  mobility  lies 
in  the  choice  of  proper  electron  effective  mass  in  eqs   (2.7)  and  (2.9)  when  energy  band 
structure  is  nonspherical .     For  silicon,  the  anisotropic  scattering  effect  comes  from  elec- 
tron mass  anisotropy   (e.g.,  m^  =  0.192  m^  and  m^  =  0.98  m^  [11])  due  to  the  ellipsoidal  con- 
duction band  structure.     Long  [14]  has  made  an  extensive  study  of  the  validity  of  Brooks- 
Herring  formula  for  n-type  silicon  for  dopant  density  less  than  10^^  cm~3  and  temperatures 
below  100  K.     He  concluded  that  if  the  electron  effective  mass  in  eq   (2.9)  was  used  as  an 
adjustable  parameter,  then  good  agreement  between  theory  and  experiment  can  be  obtained.  In 
the  present  study  [12],  we  have  also  made  a  detailed  comparison  among  the  theoretical  models 
developed  by  Brooks-Herring  [4,7],  Conwell-Weisskopf  [13],  and  Samoilovich  et  at.    [9],  and 
have  found  that  the  mobility  values  calculated  from  these  models  for  n-type  silicon  are  gen- 
erally too  large  for  10^3  <  N^^  <  10^^  cm~3  and  100  <  T  <  500  K.''"    To  correct  this  discrep- 
ancy, we  have  made  use  of  the  results  of  Long  [14],  and  Norton  et  al.    [2],  and  obtained  a 
mobility  expression,  similar  to  that  of  eq  (2.9),  for  ionized  impurity  scattering: 

Vj^  =  7.3  X  10^^  t3/2/n^  ^^^i^-  (2.10) 

Here  G(bj^)  is  identical  to  eq  (2.6)  with  the  exception  that  b  is  replaced  by  b  ,  and  m^^^  = 
0.98  m^  IS  used  in  eq  (2.7)  to  obtain  b^.     The  numerical  coefficient  in  eq  (2.10)  was  ob- 
tained directly  from  reference  [14]  for  n-type  silicon  when  the  anisotropic  scattering  ef- 
fect is  included.     The  reader  is  referred  to  the  original  paper  [14]  for  a  complete  descrip- 
tion. 

When  electron-electron  scattering  effect  is  included  in  eq  (2.10)  for  N^^  >  2  x  10^^  cm-3j 
the  formula  appears  to  yield  the  best  theoretical  fit  to  the  measured  mobility  in  the  range 
of  dopant  density  and  temperature  studied  here. 


2.3.     Neutral  Impurity  Scattering 


Because  of  the  analogy  between  a  neutral  donor  atom  and  a  hydrogen  atom,  the  mobility  for 
scattering  by  the  neutral  donors  can  be  obtained  by  proper  modification  of  the  results  for 
scattering  of  slow  electrons  by  hydrogen  atoms.     Erginsoy  first  predicted  a  temperature- 
independent  mobility  given  by  [15]: 


2TT3q3in* 
n 


N  s 


10-2  (2.11) 


This  equation  fails  to  predict  the  temperature  dependence  of  neutral  impurity  scattering  mo- 
bility observed  for  n-type  silicon  at  low  temperatures   [2].     Sclar  [16]  sought  to  improve 
eq  (2.11)  by  including  the  possibility  of  bound  states  in  the  electron-hydrogenic  impurity 
scattering  problem.     He  obtained  the  following  expression  for  the  neutral  impurity  scatter- 
ing mobility  [16] 

where 

E^,  =  1.136  X  10-15   (m'/m  )   (e  /e  )2  (2.13) 
N  n    o        o  s 


2     kT  ^/2  E  1/2 

±  (hi.)         +  ±  (-J1) 

N 


(2.12) 


tFor  details  see  Appendix  A. 


5 


10 


=  1111 


CO 


10 


E 
o 


5  10 


CD 


5  10  - 


o 

00 


lOV- 


1  10 


10 


nTTTT]  1  I  I  I  I  I  ||lll|  TTT 


I  I  1 1 


J_L 


I  I  IIIIIU 


\ 

W 


T  =  300  K 

  SCLAR 

  ERGINSOY 


\ 
\  \ 


W 


w 


W 


Hi 


10 


14 


J5 


.16 


J7 


10  10  10  10 

DONOR  DENSITY  (cm-^) 


18 


10 


19 


Figure  2.  Neutral  impurity  scattering  mobility  as  a  function  of  donor  density  computed 
from  eq  (2.11)    (due  to  Enginsoy)  -  eq  (2,12)   (due  to  Sclar) . 
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and  )i,„  is  given  by  eq  (2.11). 

Nil 

Equation  (2.12)  predicts  that        varies  as  T^^^  for  kT  above  the  binding  energy,  E^.  The 
experimental  data  given  by  Norton  et  al.    [2]  for  n-type  silicon  show  such  a  dependence  up  to 
50  K.     Our  mobility  and  resistivity  calculations  using  eqs  (2.11)  and  (2.12)  show  that  the 
latter  produces  a  better  fit  with  the  experimental  data  than  the  former.     For  comparison, 
the  neutral  impurity  scattering  mobilities  computed  from  eqs  (2.11)  and  (2.12)  are  shown  in 
figure  2  for  T  =  300  K.     Note  that  the  mobility  values  predicted  by  eq  (2.11)  are  somewhat 
smaller  than  those  of  eq  (2.12). 

2.4.     Effect  of  Electron-Electron  Scattering 

The  mobility  formula  given  in  eq  (2.10)  neglects  the  effect  of  electron-electron  (e-e)  scat- 
tering on  the  ionized  impurity  scattering  mobility.  Although  e-e  scattering  does  not  affect 
the  current  density  directly  since  it  cannot  alter  the  total  momentum,  it  tends  to  randomize 
the  way  in  which  this  total  momentum  is  distributed  among  electrons  with  different  energy. 
When  the  scattering  mechanism  is  such  as  to  lead  to  a  nonuniform  distribution,  e-e  scatter- 
ing gives  rise  to  a  net  transfer  of  momentum  from  electrons  which  dissipate  momentum  less  ef 
ficiently  to  those  which  dissipate  more  efficiently,  resulting  in  an  overall  greater  rate  of 
momentum  transfer  and  lower  mobility. 

On  the  basis  of  the  above  argument,  it  is  obvious  that  the  size  of  the  effect  of  e-e  scatter 
ing  on  the  mobility  is  a  function  of  the  energy  dependence  of  the  relaxation  time.  Thus, 
for  neutral  impurity  scattering  where  the  relaxation  time  is  independent  of  energy,  the  mo- 
bility is  not  affected  by  e-e  scattering.     Ionized  impurity  scattering  would  be  expected  to 
be  much  more  affected  than  lattice  scattering  since  in  the  former  case        is  proportional  to 
E^'^,  while  in  the  latter  it  is  proportional  to  E"-^'^. 

Luong  and  Shaw  [10]  have  analyzed  the  effect  of  e-e  scattering  on  the  ionized  impurity  mobil 
ity  using  a  single-particle-like  approximation  from  the  time-independent  Hartree-Fock  theory 
They  have  shown  that  with  the  correction  for  e-e  scattering,  the  Brooks-Herring  formula  is 
reduced  by  a  factor  which  can  be  expressed  in  closed  form  as  [10] 

m'^  =  (1  -  e'^^h)  (2.14) 

where  N^  is  the  ionized  impurity  density,  n  is  the  electron  density,  and  is  given  by  eq 
(2.10).  For  uncompensated  n-type  silicon,  the  density  of  ionized  donor  impurities,  N^,  is 
equal  to  the  density  of  conduction  electrons,  n.     Thus,  eq  (2.14)  reduces  to 

=  (1  -  e-1)  y^  =  0.632  y^  (2.15) 

The  factor  in  eq  (2.15)  is  in  good  agreement  with  a  previous  prediction  [17]  based  on  the 
Boltzmann  theory.     For  n-type  silicon,   the  effect  of  e-e  scattering  on  the  ionized  impurity 
scattering  becomes  important  for  N^^  >  2  x  10^^  cm~^,  as  will  be  discussed  later. 

The  effect  of  e-e  scattering  on  the  lattice  mobility  has  also  been  discussed  in  several 
classical  papers  [4,18,19].     It  can  be  shown  [18,19]  that  e-e  scattering  reduces  the  lattice 
scattering  mobility  by  a  maximum  of  12  percent.     This  factor  will  be  used  in  the  present  cal 
culations . 

3.     CALCULATIONS  OF  ELECTRON  MOBILITY  AND  RESISTIVITY  IN  iV-TYPE  SILICON 

3.1.     Electron  Mobility  Calculations 

We  now  discuss  the  mobility  calculations  for  n-type  silicon  in  the  range  of  dopant  density 
from  10^3  to  10^^  cm"^  and  temperature  from  100  to  500  K. 

The  combined  mobility  due  to  both  lattice  and  ionized  impurity  scattering  contributions  can 
be  calculated  according  to  the  mixed-scattering  formula  [18]: 


y^^  =  y^  [1  +  x2  (CiX  cosX  +  sinX  (SiX  -  |-)  ] 


(3.1) 


where        is  the  lattice  mobility, 

X2  =  6y^/y^, 


(3.2) 


and        is  given  by  eq  (2.10). 

CiX  and  SiX  are  the  cosine  and  sine  integrals  of  X,  respectively.     Note  that  eqs  (3.1)  and 
(3.2)  are  applicable  for  donor  densities  less  than  2  x  10^^  cm~^  where  the  effect  of  e-e 
scattering  is  negligible.     For  dopant  densities  greater  than  2  x  10^^  cm~^,  the  effect  of 
e-e  scattering  is  incorporated  empirically  in  eq  (3.2)  as  follows: 

(i)     High  dopant  density  range  (2  x  10^     <  N^^  <  10^^  cm"^)  . 

In  this  dopant  density  range,  experimental  evidence  indicates  that  the  full  effect  of  e-e 
scattering  should  be  taken  into  account  in  the  mobility  calculations,  and  eq  (3.2)  is  re- 
placed by: 

X"2  =  6  ^s'^/v'^  ,  (3.3) 

where 

=  0.88  y^,  (3.4) 


and , 


y^'  =  0.632  y^.  (3.5) 


(ii)     Intermediate  dopant  density  range  (2  x  10^^  £       £  ^  x  10^''  cm  ^)  . 

In  this  dopant  density  range,  the  effect  of  e-e  scattering  on  both  lattice  and  ionized  impu- 
rity scattering  mobilities  is  increased  gradually  with  increasing  dopant  density.  To  obtain 
the  best  fit  to  the  measured  mobility  at  300  K  in  this  transition  range,  the  mobility  reduc- 
tion factors  R(Nq)  and  S(Nj^),  for  both  and  y^  were  derived  empirically,  assuming  a  linear 
dependence  on  the  donor  density.     The  results  are  given  by: 

Vl  =  R(N^)  y^^,  (3.6) 


where 


and 


where 


R(Nj^)  =  1.013  -  6.62  X  10-19  x  (3.7) 


\s[  =  S(Nj^)  yj,  (3.8) 


S(Nj^)  =  1.04  -  2.04  X  10-18  X         .  (3.9) 

The  mobility  ratio  in  eq  (3.2)  is  now  replaced  by 

X'2  =  6  y^/pi  (3.10) 

Note  that  the  mobility  reduction  factors,  R(Nj,)  and  S(Np),  in  eqs  (3.7)  and  (3.9)  decrease 

linearly  from  1  to  0.88  and  0.632,  respectively,  as  the  dopant  density,  N  ,  increases  from 
2  X  10 16  to  2  X  10 1 '7  cm"3. 

Equations  (3.1)   through  (3.10)  allow  the  calculations  of  y      for  lO^^  1  10^^  cm-^,  and 

100  <  T  <  500  K.  When  neutral  impurity  scattering  is  included,  the  total  electron  mobility 
may  be  computed  from  the  expression 


where  y^^^  is  given  by  eq  (3.1)  and  V-^  is  the  mobility  due  to  neutral  impurity  scattering  as 
given  in  eq  (2.12).     A  sum  of  the  reciprocal  mobilities  for  this  case  is  expected  to  yield  a 
considerably  better  approximation  than  it  does  for  the  reciprocal  sum  of        and  y  because 
neutral  impurity  scattering,  being  energy  independent,  does  not  affect  any  contributions 
from        and  y^. 

3.2.     Ionization  of  the  Donor  Impurity 

In  this  section,  formulations  are  given  for  computing  the  ionized  donor  density  (or  electron 
density)  as  a  function  of  the  total  donor  density  in  uncompensated  n-type  silicon.     In  order 
to  compute  the  electron  mobility  and  resistivity,  it  is  necessary  to  know  the  amounts  of 
ionized  and  neutral  impurity  atoms  so  that  their  scattering  contributions  can  be  calculated. 

The  ionized  donor  densities  for  n-type  silicon  were  computed  by  solving  the  charge  balance 
equation  for  the  Fermi  energy  by  an  iteration  procedure.     Since  the  minority  carrier  density 
(i.e.,  hole  density  in  n-type  silicon)  is  negligibly  small  for  this  case,   the  charge  balance 
equation  is  simply 


N+  (3.12) 


where 


+  _ 


1+2  exp[(Ep  -        +  E^)/kT] 


(3.13) 


and  the  electron  density,  n,  is  given  by  [20] 


n  =  N^/[exp(E^  -  Ep)/kT  +  0.27],  for        <  1.3  kT  (3.14) 


where       =  2 


2tt  m*  kT 
a 


3/2 


is  the  effective  density  of  conduction  band  states   (N^  = 


3.22  X  10^^  cm  ^  for  n-type  silicon  at  300  K) .     The  temperature  dependence  of  the  density  of 
states  effective  mass  was  taken  into  account  in  accordance  with  the  results  of  Barber  [11]. 

Experimental  evidence  exists  which  shows  that  the  donor  ionization  energy,  e^^,  is  not  con- 
stant, but  decreases  with  increasing  dopant  density.     Hall  coefficient  measurements  by  Pear- 
son and  Bardeen  [21]  and  more  recently  by  Penin  et  al.    [22]  in  heavily-doped  silicon  from  4 
to  300  K  show  no  evidence  of  an  ionization  energy  at  impurity  densities  greater  than 
3  X  10^^  cm~^.     For  n-type  silicon  doped  with  phosphorus  impurities,  the  dependence  of  donor 
ionization  energy  on  the  dopant  density  is  [18] 

=  e^(o)  -  a  N^^^/s  (3.15) 

where  (^^(o)  =  0.045  eV  is  the  ionization  energy  of  the  phosphorus  donor;  with  a  =  3.1  x  10~^, 
eq  (3.15)  gives  a  zero  ionization  energy  for  N^^  =  3  x  10^^  cm~^. 

3.3.     Resistivity  Analysis 

Measurements  of  resistivity  in  silicon  have  been  reported  by  previous  investigators  [23-35]. 
Irvin  [25]  first  showed  complete  resistivity  versus  dopant  density  curves  for  both  n-  and  p- 
type  silicon  at  300  K,  using  mostly  previously  published  data.     Recently  Mousty  et  al.  [26] 
reported  the  resistivity  versus  phosphorus  density  in  n-type  silicon.     The  electrical  prop- 
erties of  heavily-doped  silicon  were  also  reported  by  Chapman  et  at.   [28].     Most  of  these  ef- 
forts focused  on  resistivity  measurements  near  room  temperature.     A  theoretical  analysis  of 
the  resistivity  as  a  function  of  dopant  density  for  n-type  silicon  has  been  reported  by  Rode 
[3]  and  Norton  et  al.    [2];  their  calculations  of  resistivity  were  also  confined  to  room  tem- 
perature.    In  this  work,  we  have  extended  the  theoretical  calculations  of  resistivity  as  a 
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function  of  dopant  density  for  n-type  silicon  for  10^^  £  ^  10  cm  ^,  and  100  £  T  <_  500  K. 
The  resistivity  is  calculated  from  the  expression 

P  =  (3.16) 

If  both  p  and  n  are  measured  independently,  then  the  electron  mobility  may  be  deduced  from 
eq  (3.16).     On  the  other  hand,  if  y    and  n  are  computed  from  the  theoretical  expressions  dis- 
cussed above,  it  is  then  possible  to  compare  between  theoretically  computed  and  experimental- 
ly measured  resistivity.     We  shall  discuss  this  further  in  section  4.2. 

4.     RESULTS  AND  DISCUSSIONS 

In  this  section,  calculations  of  partial  mobilities,  electron  mobility,  electron  density, 
and  resistivity  as  functions  of  dopant  density  and  temperature,  based  on  the  theoretical  for- 
mulas discussed  in  the  preceding  section,  are  described.     Results  of  the  resistivity  measure- 
ments and  electron  mobility  data  for  phosphorus -doped  silicon  are  presented. 

4.1.     Electron  Mobility  vs.  Dopant  Density  and  Temperature 

The  mobility  formulas  used  in  the  present  calculations  are  listed  in  table  1.     As  noted  in 
this  table  the  mobility  calculations  carried  out  in  three  distinct  dopant  density  ranges  are 
based  on  our  analysis  of  the  electron  mobility  data  versus  dopant  density  at  300  K  and  the 
resistivity  data  for  phosphorus-doped  silicon  taken  at  temperatures  between  100  and  500  K. 
The  mobility  formulas  in  column  (I)  of  table  1  may  be  referred  to  as  the  traditional  mobil- 
ity formulations  which  neglect  the  effect  of  electron-felectron  scattering.     The  mobility  cal- 
culations from  these  formulas  for  n-type  silicon  are  in  good  agreement  (within  ±5  percent) 
with  experimental  data  for  dopant  densities  up  to  2  x  10-^^  cm~^,  but  for  greater  dopant  den- 
sities the  calculated  mobility  is  larger  than  the  measured  mobility.     By  including  the  ef- 
fect of  e-e  scattering  in  the  mobility  formulas  shown  in  columns  (II)  and  (III) ,  we  were 
able  to  bring  into  accord  (within  ±5  percent)  the  theoretical  calculations  and  the  experi- 
mental data  at  300  K  for  dopant  densities  between  10^^  and  10^^  cm~^. 

Figure  3  shows  the  electron  mobility  as  a  function  of  donor  density  computed  at  300  K,  along 
with  experimental  data  for  phosphorus -doped  silicon.     The  solid  dots  are  taken  from  the  mo- 
bility data  of  Mousty  et  al.    [26],  and  Baccarani  and  Ostoja  [29],  corrected  to  300  K,  and 
the  open  circles  are  data  obtained  by  Buehler  et  al.    [30,31]  on  processed  phosphorus-doped 
silicon  slices.     A  planar  four -probe  structure  was  used  for  resistivity  measurements  [30] 
and  a  diffused  diode  near  the  four -probe  test  structures  was  used  for  junction  C-V  measure- 
ments of  donor  density  [31].     The  theoretical  calculations  (curve  3)  are  within  ±5  percent 
of  the  experimental  data  for  dopant  densities  from  lOl"*  to  10^^  cm  ^.     For  comparison,  the 
Irvin  mobility  curve  (curve  2)  which  was  calculated  from  the  Caughey  and  Thomas   [32]  empiri- 
cal formula  for  n-type  silicon  is  also  shown  by  a  dashed  line  in  figure  3.     To  illustrate 
the  important  effect  of  e-e  scattering,  we  have  also  shown  in  figure  3  the  theoretical  cal- 
culations of  electron  mobility  as  a  function  of  dopant  density  (see  curves  1  and  4)  with  and 
without  inclusion  of  e-e  scattering  effect  over  the  entire  range  of  dopant  density  from  10^^ 
to  10^^  cra~^.     It  is  clearly  demonstrated  that  the  mobility  values  predicted  in  curve  1  are 
smaller  than  the  measured  mobility  for  N^^  <_  2  x  10^ cm~^  when  e-e  scattering  is  included  in 
this  dopant  density  range.     On  the  other  hand,  the  mobility  values  predicted  in  curve  4  are 
larger  than  the  measured  mobility  for  Nq  >  2  x  10^^  cm~^  when  e-e  scattering  is  neglected. 
Thus,  a  transition  region  (i.e.,  2  x  1016  <^  N^^  <_  2  x  10-^^  cm~^)  exists  in  which  the  effect 
of  e-e  scattering  increases  gradually  with  increasing  dopant  density.     The  mobility  formulas 
in  table  1  were  also  applied  to  compute  the  electron  mobility  as  a  function  of  dopant  den- 
sity for  temperatures  between  100  and  500  K.     The  results  are  displayed  in  figures  4  and  5. 
Mobility  was  not  determined  directly  at  temperatures  other  than  300  K  so  no  attempt  has  been 
made  to  compare  theory  with  experiment. 

4.2.     Resistivity  vs.  Dopant  Density  and  Temperature 

Electrical  resistivity  as  a  function  of  dopant  density  for  phosphorus-doped  silicon  was  cal- 
culated using  eq  (3.16).     Values  of  \i    were  obtained  from  the  results  of  section  4.1.,  and 
the  electron  density  was  computed  by  solving  eqs   (3.12)   through  (3.14)  using  the  iteration 
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Figure  3.     Electron  mobility  versus  total  donor  density  for  n-type  silicon  at 
300  K.     (Curve  1  is  the  theoretical  calculation  with  e-e  scattering  included  for 
10^3  <        <  10^5  cm  ^,  curve  2  is  the  Irvin  curve  calculated  from  the  Caughey 
and  Thomas  empirical  formula  [30],  curve  3  is  the  exact  theoretical  calculation 
using  mobility  formulas  in  Table  1,  curve  4  is  the  theoretical  calculation  vd.th 
e-e  scattering  neglected.    Open  circles  are  NBS  data  [27],  and  solid  dots  are 
from  Mousty  et  at.    [24]  for  phosphorus-doped  silicon.) 
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Figure  4.  Electron  mobility  of  n-type  silicon  versus  temperature  for  dopant  densities 
from  10^3  to  5  x  10^^  cm"^. 
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I  procedure  mentioned  in  section  3.2.     The  results  are  displayed  in  figures  6  through  10. 

In  figure  6,  the  theoretical  curves  are  shown  for  the  percent  ionization  of  phosphorus  atoms 
versus  total  phosphorus  density  for  temperatures  between  100  and  500  K.     The  electron  densi- 
ty was  obtained  from  eq  (3.12),  which  is  a  valid  approximation  for  the  uncompensated  case. 
Figure  7  shows  the  resistivity  versus  donor  density  for  n-type  silicon  at  300  K;  the  solid 
line  is  our  theoretical  calculation,  and  the  dashed  line  is  the  Irvin  curve  [25].     The  cal- 
culated resistivity  is  5  to  10  percent  smaller  than  the  resistivity  given  by  the  Irvin  curve 
I   for  the  same  value  of  dopant  density;  this  is  consistent  with  recently  published  resistivity 
data  [26,31]. 

In  order  to  compare  theory  with  experiment  for  temperature  other  than  300  K,  resistivity  mea 
surements  were  made  between  100  and  500  K  on  seven  phosphorus-doped  silicon  wafers  with  den-r 
sities  from  1.2  x  lO^^*  to  2.5  x  10^^  cm"^. 

The  resistivity  measurements  were  performed  using  a  planar  four -probe  test  structure  which 
was  designed  for  bulk  resistivity  measurements  and  was  fabricated  using  the  bipolar  transis- 
tor process.     The  test  cells  used  for  the  resistivity  measurements  were  obtained  from  the 
same  silicon  wafers  fabricated  for  the  mobility  measurements  at  300  K  (see  figure  3) .  For 
temperature  dependent  measurements,  the  test  cell  containing  the  four -probe  structure  was 
mounted  on  a  TO-5  header,  and  a  temperature-sensing  diode  was  also  mounted  next  to  the  cell 
for  measuring  the  temperature.     The  TO-5  header  was  then  mounted  inside  a  specially-designed 
cryostat  where  temperatures  can  be  varied  from  -191  to  400°C  with  a  maximum  heating  rate  of 
7  K/s  [32]. 

Figure  8  shows  the  plot  of  resistivity  as  a  function  of  temperature  for  phosphorus-doped  sil 
icon  with  dopant  densities  from  1.2  x  10^^  to  2.5  x  10-^^  cm~^;  the  solid  lines  are  the  theo- 
retical calculations  and  the  solid  dots  are  the  measured  resistivity.  Figure  9  displays  the 
calculated  resistivity  versus  temperature  for  dopant  densities  from  5  x  10^^  to  10^^  cm~^  in 
more  regular  steps.  Figure  10  shows  the  calculated  resistivity  curves  as  a  function  of  dop- 
ant density  for  temperatures  between  100  and  500  K.  The  results  show  that  the  resistivity 
depends  strongly  on  temperature  for  N  <  10^^  cm~^  where  lattice  scattering  is  dominant. 
However,  the  resistivity  is  nearly  independent  of  temperature  for  N^^  >  10^^  cm~^  where  mixed 
scattering  prevails. 

5.  CONCLUSIONS 

An  improved  theoretical  model  has  been  formulated  for  computing  electron  mobility  in  n-type 
silicon  as  functions  of  dopant  density  and  temperature.     In  addition  to  considering  the  con- 
tributions from  scattering  by  lattice  phonons,  ionized  impurities,  and  neutral  impurities, 
this  model  takes  into  account  the  effect  of  electron-electron  scattering  on  both  lattice  and 
ionized  impurity  scattering  mobilities  for  dopant  densities  greater  than  2  x  10^^  cm~^. 

From  this  study,  it  was  found  that  the  influence  of  e-e  scattering  is  negligible  for  <_ 

2  X  10^^  cm~^,  but  is  significant  for  N    >^  10-^^  cm~^.     To  allow  a  smooth  change  in  the  cal- 
culated mobility  values  across  the  boundaries  of  the  transition  region  (i.e.,  2  x  10^^  —  '^D 
<_  2  X  10-^ cm~^),   the  effect  of  e-e  scattering  is  assumed  to  increase  linearly  with  dopant 
density  in  this  region.     The  theoretical  predictions,  based  on  such  an  assumption,  are  with- 
in ±5  percent  of  the  measured  mobilities  for  10-^^  £       £  10^^  cm~^,  and  T  =  300  K. 

Resistivity  analysis  for  phosphorus-doped  silicon  shows  that  for  dopant  densities  less  than 

3  X  10^^  cm~^,  the  calculated  resistivity  values  are  within  ±7  percent  of  the  experimental 
data  for  temperatures  between  100  and  500  K. 
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Figure  8.     Resistivity  versus  temperature  for  seven  phosphorus-doped  silicon 
slices.     (Solid  lines  are  the  theoretical  calculations  and  dots  are  the 
experimental  data.) 
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Figure  9.  Theoretical  calculations  of  resistivity  versus  temperature  for  n-type 
silicon  with  donor  densities  from  5  x  10^^  to  10^^  cm~^ . 
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Figure  10.  Theoretical  calculations  of  resistivity  versus  total  donor  density 
for  n-ty-pe  silicon  with  temperature  as  a  parameter  (100<T<500  K) . 
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APPENDIX  A:     COMPARISON  OF  THE  EXISTING  MODELS  FOR  IONIZED  IMPURITY  SCATTERING 


In  addition  to  the  Brooks-Herring  (BH)  mobility  formula  already  discussed  in  section  2.2.  of 
the  main  text,  we  now  proceed  to  discuss  other  existing  models  for  computing  the  ionized  im- 
purity scattering  mobility  in  n-type  silicon.     These  include  the  Conwell-Weisskopf  (CW)  for- 
mula [13]  and  the  formulas  derived  more  recently  by  Samoilovich,  Korenblit,  and  Dakhovskii 
(SKD)  [9]. 

(i)     Conwell-Weisskopf  (CW)  Formula: 


Using  the  Born  approximation  and  assuming  that  the  impurity  potential  is  a  bare  Coulomb  po- 
tential (i.e.,  V(r)  =  — ^) ,  Conwell  and  Weisskopf  derived  a  mobility  expression  which  reads 
[13]: 


e  r 
s 


I(CW) 


128(2^)1/2  ^  2^3/2 
s 

qVl/2 


In  [1  +  ( 


127r  e  kT  2 

'     ^  Nj-2/3 


-) 


-1 


4.55  X  10 


17  t3/2 


,  ,   4.42  X  IQS  t2^ 

In   [1  +   ;  ] 

Nj2/3 


-1 


(A.l) 


The  numerical  values  in  the  brackets  of  eq  (A.l)  are  obtained  by  substituting  m    =  0.26  m 
and  e    =  11.7  x  8.852  x  10~^^  F/cm  in  the  above  equation.     The  theoretical  values  of  cal- 
culated from  eq  (A.l)  are  shown  in  figure  A.l  (see  curve  4). 

(ii)     The  SKD  Formula: 

Samoilovich  et  al.  [9]  have  derived  a  theoretical  model  for  computing  ionized  impurity  scat- 
tering mobilities  in  n-type  silicon  taking  into  account  the  effect  of  scattering  anisotropy. 
The  results  of  their  mobility  formulations  can  be  summarized  as  follows. 


The  average  relaxation  times  along  the  transverse  and  longitudinal  axis  of  the  ellipsoidal 
conduction  band  valley  of  silicon  are  given,  respectively,  by  [2]: 


<T^>  =  5.48  X  iQll  (g^  +  1)/{N^  X  [In  (b^)  -  2.34  +  7.88/bj^]} 


(A. 2) 


and , 


<T^>  =  1.98  X  10I2  (g^  +  1)/{N^  X  [In(b^)  -  1.75  +  4.64/b^]} 


(A.  3) 


where 


g^  =  0.03  X  [logio  (h^)  -  1] 


g^  =  0.192  -  0.335/logio  (b^) 
b^  =  1.52  X  10^5  x2/n' 


n'  =  n  +  n-(l  -  n/N^) ,  for        =  0 


(A. 4) 
(A. 5) 
(A. 6) 
(A. 7) 


The  mobility  is  then  computed  from  the  following  expression: 


^^KSKD)  =  f  ^  ^)  =  ^  ''t>  -"'^^ 


(A.  8) 


Here  <t  >  and  <t  >  are  the  Maxwellian  average  relaxation  time  components  given  by  eqs  (A. 2) 
and  (A. 5),  respectively.     The  mobility  versus  dopant  density,  as  calculated  from  eq  (A. 8), 
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Figure  A-1.     Ionized  impurity  scattering  mobility  versus  total  donor  density  as 
computed  from  different  theoretical  models  at  300  K.     (Curve  1  is  calculated  from 
eq  (2.10),  curve  2  from  eq  (2,9),  curve  3  from  eq  (A-8) ,  and  curve  4  from  eq  (A-1).) 
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is  shown  in  figure  A.l  (i.e.,  curve  3).     For  comparison,  the  mobility  curves  calculated  from 
the  Brooks-Herring  formula  (i.e.,  eq  (2.9))  and  the  modified  formula  (i.e.,  eq  (2.10))  dis- 
cussed in  the  main  text  are  also  included  in  figure  A.l  (see  curves  2  and  1).     Note  that  cor- 
rections due  to  e-e  scattering  for  N^^  >  2  x  10^^  cm~^  have  also  been  taken  into  account  for 
all  four  curves  shown  in  figure  A.l.     The  results  show  that  for  N^^  <  10^^  cm~^  the  mobility 
values  predicted  by  the  three  existing  models  are  in  reasonable  agreement  (within  ±20  per- 
cent), but  are  25  to  50  percent  larger  than  those  predicted  by  the  modified  formula  (i.e., 
eq  (2.10)).     However,  for  N^^  >  10^*^  cm~^,  the  mobility  values  calculated  by  both  the  SKD  and 
BH  formulas  are  considerably  larger  (about  50  to  100  percent)   than  those  from  eq  (2.11). 

In  conclusion,  the  present  results  suggest  that  the  mobility  formulations  for  ionized  impu- 
rity scattering  derived  by  Conwell-Weisskopf ,  Brooks-Herring,  and  Samoilovich  et  at.  appear 
to  overestimate  the  mobility  for  n-type  silicon,  particularly  for  N^^  >  10^^  cm~^.     One  possi- 
ble reason  for  this  discrepancy  may  be  attributed  to  the  anisotropic  scattering  effect  re- 
sulting from  electron  mass  anisotropy  in  n-type  silicon,  as  pointed  out  in  the  main  text. 
For  p-type  silicon  where  hole  effective  masses  are  nearly  isotropic,  the  above  mentioned  dif- 
ficulty would  not  be  expected  to  occur. 
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APPENDIX  B:     DIELECTRIC  SCREENING  LENGTH  AND  DEGENERACY  TEMPERATURE  VS.  DOPANT  DENSITY 


In  this  appendix,  we  calculate  two  important  physical  parameters  which  are  relevant  to  the 
mobility  calculations  discussed  in  the  main  text  and  Appendix  A.     These  two  parameters  are 
the  dielectric  screening  length,  R  ,  and  the  degeneracy  temperature,  T , .      For  the  sake  of 
simplicity,  we  consider  here  the  case  of  phosphorus-doped  silicon  in  which  one  extra  elec- 
tron is  provided  by  each  phosphorus  donor  atom  which  then  becomes  positively  charged  when 
ionized.     Let  n  be  the  density  of  electrons  in  the  conduction  band,  and  V(r)  the  electro- 
static potential.     In  the  presence  of  the  self-consistent  field  due  to  V(r),  the  ordinary 
space-independent  Fermi  energy,  E  ,  must  be  replaced  by  E    +  qV(r),  and  V(r)  can  be  obtained 
by  solving  Poisson's  equation.     The  result  is  given  by  [71: 

V(r)  =  C-^)  exp  (-r/R^)  (B.l) 


s 

where 

R    =  O   (7^)  (B.2) 
n        q  47rn 

For  n-type  silicon,  eq  (B.2)  reduces  to 

R    =  23.6  t1/2  (cm)  (B.3) 

n 


which  shows  that  the  dielectric  screening  length  is  inversely  proportional  to  the  square 
root  of  the  electron  density.     Equation  (B.3)  was  used  to  calculate  the  dielectric  screening 
length  versus  dopant  density  for  n-type  silicon  for  temperatures  between  100  and  500  K.  The 
results  are  shown  in  figure  B.l.     For  comparison,  the  interatomic  spacing  of  the  silicon  lat- 
tice is  also  included  in  this  figure.     Note  that  the  dielectric  screening  of  conduction  elec- 
trons may  be  considered  "complete"  if  R    >>  d  . ,  where  d  .  is  the  spacing  of  silicon  atoms. 
For  n-type  silicon,  this  condition  is  satisfilS  for  N    <  ^  x  10^^  cm"^  and  T  >  100  K.  Thus, 
calculations  of  ionized  impurity  scattering  mobility  By  the  Brooks-Herring  formula,  as  dis- 
cussed in  the  main  text,  should  be  valid  for  N^^  <^  3  x  10 cm~^  and  T  >  100  K.     However,  for 
higher  temperatures  (e.g.,  T  >_  250  K) ,  the  validity  of  the  Brooks-Herring  formula  may  be  ex- 
tended to  N    '\'  10 1^  cm~^,  since  the  dielectric  screening  length  also  increases  with  increas- 
ing temperature. 

The  degeneracy  temperature  of  the  electron  system  in  n-type  silicon  is  calculated  from  the 
following  expression  [7]: 

h2        3n  2/3 

T,  =  (-V)   (^)  (B.4) 

8m\ 
n 


Equation  (B.4)  was  used  to  compute  the  degeneracy  temperature  as  a  function  of  electron  den- 
sity for  n-type  silicon,  and  the  result  is  shown  in  figure  B.2.     The  curve  may  serve  as  a 
demarcation  line  which  separates  the  degenerate  region  (i.e.,  above  the  curve)  from  the  non- 
degenerate  region  (i.e.,  below  the  curve)  as  shown  in  the  figure.     This  plot  is  useful  for 
identifying  the  validity  of  the  statistics  used  in  the  mobility  calculations.     For  example, 
at  T  =  100  K,  nondegenerate  statistics  may  be  used  in  the  mobility  calculations  for  dopant 
densities  up  to  3  x  10^^  cm"^  and  this  may  be  extended  to  7  x  lO-^^  cm"^  for  T  =  300  K. 
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Nolcm-^) 


Figure  B-1.  Dielectric  screening  length  versus  donor  density  for  100<T<500  K  for 
phosphorus-doped  silicon. 
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Figure  B-2.  Degeneracy  temperature  as  a  function  of  electron  density  for  n-type 
silicon. 
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